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Effects of electron-electron and spin-orbit interactions on the ground-state magnetic configuration 
and on the corresponding thermoelectric and spin thermoelectric properties in zigzag nanoribbons of 
two-dimensional hexagonal crystals are analyzed theoretically. Thermoelectric properties of quasi¬ 
stable magnetic states are also considered. Of particular interest is the influence of Coulomb and 
spin-orbit interactions on the topological edge states and on the transition between the topological 
insulator and conventional gap insulator states. It is shown that the interplay of both interactions 
has also a significant impact on transport and thermoelectric characteristics of the nanoribbons. 
The spin-orbit interaction additionally determines an in-plane magnetic easy axis. Thermoelectric 
properties of the nanoribbons with in-plane magnetic moments are compared with those of the 
nanoribbons with edge magnetic moments oriented perpendicularly to their plane. Nanoribbons 
with ferromagnetic alignment of the edge moments are shown to reveal spin thermoelectricity, in 
addition to the conventional one. 
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I. INTRODUCTION 

There is currently a great interest in strictly two- 
dimensional (2D) crystalline materials, like graphene^i 
silicene)^!^ germanene^ and others. A common feature 
of these materials is a 2D hexagonal lattice of atoms - 
carbon, silicon, and germanium, respectively. However, 
both silicene and germanene have buckled atomic struc¬ 
ture, where the two triangular sublattices are slightly dis¬ 
placed along the normal to the atomic plane, whereas 
all carbon atoms in graphene are localized in a com¬ 
mon plane. This buckling of atomic structure leads to 
new properties of silicene and germanene,—^— which make 
them different from graphene. 

Electronic states around the charge neutrality points 
of these materials are usually described by the relativis¬ 
tic Dirac modeli^ii^ This model, however, neglects the 
electron-electron interaction. Moreover, some deviations 
of the exact electronic structure from the conical disper¬ 
sion relation characteristic of the Dirac model appear for 
states distant from the Dirac points. In other words, the 
exact electronic structure is more complex, and either 
tight binding models or ab-inito techniques are usually 
employed to calculate the corresponding band structure, 
especially in the case of nanoribbonsIt is already 
well known that the Coulomb interaction between elec¬ 
trons leads to edge magnetism in nanoribbons with zigzag 
edgesii^ In the ground state of the system, magnetic mo¬ 
ments localized at one edge are opposite to those at the 
other edgeji2r2i In some situations, however, one can sta¬ 
bilize a quasi-stable ferromagnetic state, where all edge 
moments are oriented in paralleli^^— Spin-orbit interac¬ 
tion gives rise to a magnetic anisotropy, and when this 
coupling is sufficiently strong, an in-plane orientation of 
the edge moments is energetically more favorable<^ Thus, 


the interplay of spin-orbit and Coulomb interactions is of 
fundamental importance for the edge magnetism. 

The spin-orbit interaction also leads to topologically 
protected edge states in the spin-orbit gap?^ and these 
states are responsible for topological-insulator properties 
of these materialsj ^^i^^ However, Coulomb interaction has 
a significant impact on these properties. Accordingly, the 
interplay of spin-orbit and Coulomb interactions is im¬ 
portant not only for the edge magnetism, but also for 
the topological state of these materials. Therefore, in 
this paper we consider this issue in the case of zigzag 
nanoribbons of two-dimensional buckled hexagonal crys¬ 
tals, and especially the role of Coulomb interaction and 
spin-orbit coupling. We focus on the thermoelectric and 
spin thermoelectric effects of these materials in the bal¬ 
listic transport regime, and on their dependence on the 
Coulomb and spin-orbit interactions. First, we calcu¬ 
late energy of various magnetic configurations in zigzag 
nanoribbons, and find that the antiparallel configura¬ 
tion with in-plane magnetic moments is the ground state 
structure for large enough spin-orbit coupling. For small 
spin-orbit coupling, the in-plane and perpendicular-to- 
plane configurations are degenerate. Then, we calculate 
the electronic states, with a special emphasis on the edge 
states responsible for the topological properties. Having 
determined the magnetic state and also the electronic 
edge states, we calculate thermoelectric properties with 
the main focus on the role of topological edge states and 
transition from the topological to conventional insulator 
state. 

Generally, the thermoelectric and especially the spin 
thermoelectric effects in systems of reduced dimensional¬ 
ity are currently of great interest mainly due to a 
hope to find an efficient way to convert dissipated heat 
into electrical energy. In fact, thermoelectric properties 
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of nanoribbons have been studied theoretically in a cou¬ 
ple of papersHowever, the role of Coulomb interac¬ 
tion and topological edge states in the gap, has not been 
studied thoroughly enough yet. In a recent paper— we 
have analyzed the influence of topological states on the 
thermoelectric properties, with the main focus on the role 
of a staggered exchange field^^— and of an electric field 
perpendicular to the atomic plane. Here we consider in 
more details the role of Coulomb interaction, which has 
a significant impact on the topological properties— 

In section 2 we present the Hamiltonian used to de¬ 
scribe the 2D materials under consideration, and also in¬ 
troduce briefly the method used to calculate transmission 
function and then all the thermoelectric and transport 
coefficients. In section 3 we present numerical results 
on the edge magnetism, topologically protected states, 
and thermoelectric properties of the corresponding zigzag 
nanoribbons. Summary and final conclusions are in sec¬ 
tion 4. 


II. THEORETICAL DESCRIPTION 
A. Tight-binding Hamiltonian 

The zigzag nanoribbons of hexagonal 2D crystals, like 
silicene or germanene, which are analyzed in this paper 
consist of N zigzag atomic chains. The system can be 
described by Hamiltonian of the following general form: 

H = Hth + H,o + He-e, (1) 

where the first term, Hth, is the tight-binding Hamil¬ 
tonian, the term iJgo represents the spin-orbit coupling, 
while the last term, He-e, stands for the electron-electron 
Coulomb interaction. 

The tight-binding Hamiltonian is assumed in the form 

Hth = —t ( 2 ) 

(iJ) 

where (cja) is the creation (annihilation) operator 
for an electron with spin cr =f,at the lattice point i 
(j). The spin quantization axis is taken along the z-axis, 
which is normal to the atomic plane. The tight-binding 
term is restricted here to electron hopping between 
nearest-neighbor sites, with the corresponding hopping 
parameter t. For instance, t = 1.6 eV in silicene<^ 

The intrinsic spin-orbit interaction in Eq. o follows 
from electron hopping between next-nearest-neighbors in 
the hexagonal lattice, and can be written in the form^i 

Hso = *Aso Cj^ — Cjl^ , (3) 

where Ago is the spin-orbit interaction parameter, while 
Vij = 1 (vij = — 1) when the hopping path from second- 
neighbor sites i to j in the hexagonal lattice is clockwise 
(anticlockwise) with respect to the positive z-axis. 


The electron-electron interaction will be taken into ac¬ 
count in the Hubbard form within the mean-field approx¬ 
imation (MFA). The form of the MFA Hubbard Hamil¬ 
tonian depends on whether the spin quantization axis 
is parallel to the local magnetization orientation or not. 
We assume, similarly as in Ref.@, that in a hypothetical 
state the magnetic moments are tilted from the normal 
orientation by an angle a, and will calculate the cor¬ 
responding total energy. The ground state configuration 
will be thus the one with the lowest energy. Assuming the 
quantization axis parallel to the local magnetic moments, 
the Hubbard term in the mean-field approximation can 
be written in the form 

He-e = (<4) + <4 (<t) “ (<t) (<4)) ’ ( 4 ) 

i 

where and are the particle number operators for a 
site i and for electrons with spin f and ]. with respect to 
the spin quantization axis tilted from the normal to the 
plane towards a nanoribbon edge by an angle a. In turn, 
U is the on-site Coulomb repulsion parameter, U > 0. 
We also note, that the Coulomb interaction of electrons 
at different lattice sites is omitted. 

The MFA Hubbard Hamiltonian should be now writ¬ 
ten in the quantization axis common for all terms of the 
Hamiltonian. This can be done by rotating the quanti¬ 
zation axis back to the orientation normal to the plane, 
which is achieved by the transformation^ 



where a = 1 (ft = —1) for a =t (cr =),). 

The spin-dependent mean number of electrons, 
at the site i, required to evaluate the Hubbard interaction 
term o, is determined self-consistently from the band- 
structure of the corresponding nanoribbon, 

27r/o 

j f[Eik),fj]dk, (6) 

bands q 

where /(E, /i) is the Fermi-Dirac distribution function, 
with fj, denoting the chemical potential, while nf^{k) is 
the spin-dependent electron density at site i for a Bloch- 
wave of energy E{k). The parameter a denotes the length 
of the elementary cell of a zigzag nanoribbon. The chem¬ 
ical potential of a pristine nanoribbon is equal to the 
corresponding Fermi energy Ep. By doping with donor 
or acceptor impurities, or by applying an external gate 
voltage, one can shift up or down the Fermi level, and its 
position with respect to Ep will be described by /r — Ep. 
In the following f and 4 . will denote spin projection on 
the quantization axis parallel to the orientation of the 
magnetic moments. 

To identify the edge-states, we evaluate the mean value 
of ^ = 2(y — 2 /o)/ty, where the y-axis is in the nanoribbon 
plane and normal to the nanoribbon axis, yo is the y- 
coordinate of the nanoribbon center, while w is the width 
of the nanoribbon. For ^ = — 1 (^ = 1), the states are 
localized at the left (right) edge of the nanoribbon. 
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B. Electronic transmission and general formulas 
for thermoelectric coefficients 

Ballistic electronic transmission in quasi-one- 
dimensional structures is usually determined in terms 
of the Green’s functions formalism and an expression 
derived by Caroli et ol— An alternative approach is 
based on matching wave functions in the scattering 
region to the Bloch modes of ideal bulk leads4i The 
total transmission, Ttot(i'), as a function of electron 
energy E is then defined as the sum of transmissions 
for all incident Bloch waves for a given wave-number 
k, and is a piece-wise constant function equal to the 
number of energy bands for a given energy E. In turn, 
the spin transmission, Tspin, is equal to the sum of the 
mean values of the Pauli operator dz. More detailed 
description of the transmission calculation is presented 
in Ref. [i^. The total transmission, Ttot, and spin 
transmission, Tspin, are related to the spin-dependent 
transmission, T„ (a =t, 4 .), as Ttot = + Tj, and 

Tspin = T^-Ti, respectively. 

Knowledge of the transmission function Ta{E) allows 
to calculate all transport characteristics in the linear re¬ 
sponse regime, including also the thermoelectric prop¬ 
erties. The spin resolved transport and thermoelectric 
coefficients can be then determined from the moments 
Lma of the electronic current series expansion around 
the Fermi level, 

CXD 

Lm. = \ J T,{E) {E - m)™ dE. (7) 


When the spin thermopower is absent, one can write 
the conventional Seebeck coefficient S {Sc —>■ S) in the 
standard form S = —(l/|e|T)(Ti^ -|- Ti^/Tot -I- Tqj,) = 
-(l/|e|T)(Ti/To), with T„ = T„| -|- for n = 0,1,2. 
In turn, the electronic contribution k to the thermal con¬ 
ductance is given by k = (1 /T)(T2 — L\/Lq). 

Finally, one can define the dimensionless figures of 
merit for the conventional and spin thermoelectric effects, 


ZTc = 


ZT, = 


SlGT 

K -I- Kph ’ 
K + Kph 


( 10 a) 

( 10 b) 


where Kph is the phonon contribution to the heat con¬ 
ductance. In the absence of spin thermoelectricity, only 
the dimensionless figure of merit for the conventional 
thermoelectricity is relevant {ZTc —>■ ZT), and ZT = 
52GT/(K + «:ph), 


III. NUMERICAL RESULTS 

In this section we present some numerical results on the 
basic thermoelectric properties. Before this, however, we 
consider various stable and quasi-stable magnetic states 
of the nanoribbons. Then, we consider separately the 
antiferromagnetic (AFM) states, with magnetic moments 
at one edge being antiparallel to the magnetic moments 
at the other edge), and ferromagnetic (FM) states, with 
parallel alignment of all the edge moments. 


For a piece-wise constant transmission function, calcula¬ 
tion of Lmcr can be reduced to calculation of the so-called 
incomplete Fermi-Dirac integral,— 

OO 

, I [ z'^dz 

Em{x) = — , for X > 0 , 8 ) 

ml J 1 -I- 

X 

for which there exists an efficient numerical procedure— 
The spin-dependent conductance Go- is proportional 
to the moment Too-, G^ = Lq^t where e is the electron 
charge. The total conductance G is then equal to Gy-I-Gj,, 
and the spin conductance is Gs = G^ — G^. The heat 
conductance is defined at vanishing charge current and 
can be written as k = (1/T) ~ Zia-fLoa), where 

T stands for temperature. In turn, the spin-dependent 
thermoelectric coefficient (thermopower) So- can be ex¬ 
pressed in the form 


1 Tier 

|e|T Loa 


(9) 


for a =t,4.- Then, one can define the conventional 
(charge) Seebeck coefficient (thermopower) as Sc = 
^(S”-!- -I- 54 ,) and also the corresponding spin Seebeck 
coefficient (spin thermopower) as Ss = 5 (*S't ~ 


A. Stability diagram 

We have performed numerical calculations of the total 
energy of nanoribbons with the edge magnetic moments 
tilted by an angle a from the orientation normal to the 
nanoribbon plane, for a ranging from a = 0 to a = 7 r/ 2 . 
The calculations have been carried out for both FM and 
AFM magnetic configurations. The numerical results 
clearly show that the states for a ^ 0 , 7 r /2 are neither 
stable nor quasi-stable. The only stable or quasi-stable 
states are those for a = 0 (moments normal to the plane) 
and a = 7 r /2 (moments in the nanoribbon plane). This 
applies to both FM and AFM configurations. From this 
analysis, we have constructed the stability diagram of 
the FM and AFM states. Due to the spin-orbit interac¬ 
tion, both AFM and FM configurations have lower energy 
when the edge moments are oriented in the nanoribbon 
plane than when they are along the normal to the plane, 
and the energy difference depends on the Hubbard pa¬ 
rameter JJ and spin-orbit parameter Ago. Thus, for each 
magnetic configuration one can distinguish between the 
perpendicular (T) and parallel (||) states.— In the former 
case the edge moments are normal to the nanoribbon 
plane, while in the latter one they are in the nanoribbon 
plane. As a result, one finds the following four states. 
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FM_L, AFM_L, FM|| and AFM||. Numerical results show 
that the AFM|| state is the ground state of the system, 
i.e. the state of the lowest energy for Ago > 0. 
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FIG. 1: Energy difference, AE, of the a) AFM_L, FM_L, and 
b) FM|| states and of the AFM|| gronnd state, preseirted as a 
function of the Hubbard parameter U for indicated values of 
the spin-orbit parameter Aso. Other parameters: t = 1.6 eV, 
and N = S. 


In Fig. [T] we show the energy AE of the states FM_L, 
AFM_L, and FM||, measured from the energy of the 
AFM|| ground state. The excess energy, AE, of a partic¬ 
ular state over the AFM|| ground state is presented there 
as a function of U for indicated values of the spin-orbit 
parameter Ago. One should note, that only nonmagnetic 
solutions have been found for low values of U, while mag¬ 
netic phases appear for larger values of U. ^From Fig.[TJi 
follows, that the perpendicular configuration of magnetic 
moments is not favorable in the presence of spin-orbit 
coupling, which is in agreement with the results obtained 
by Lado et ah^ Both AFM_L and FMT states correspond 
to energies higher than that of the AFM|| state, and the 
energy difference strongly increases with the spin-orbit 
parameter. As one might expect, the AMFT and AMF|| 
states are practically degenerate in the limit of small val¬ 
ues of Ago. In turn, energy of the FMT state is consid¬ 
erably higher than the energy of the AMFT state, and 
the difference increases roughly linearly with increasing 


U. Interestingly, the energy difference between the per¬ 
pendicular phases AFMT and FMT for a given U de¬ 
creases with increasing Ago, so these two configurations 
become almost degenerate for strong spin-orbit coupling 
and moderate values of the Hubbard parameter U. Re¬ 
markably different behavior can be observed for the FM|| 
state, as shown in Fig. [Hr. For high values of Ago and 
moderate values of U, this state is practically degenerate 
with the AFM|| one, but the energy of FM|| state no¬ 
tably increases when the spin-orbit interaction becomes 
small. Finally, in the limit of Ago —?► 0, the states AFM|| 
and AFMT are fully degenerate for all relevant values of 
U (Fig. (Hi), whereas the FM states are quasi-stable and 
correspond to higher values of energy. 

Below we will analyze transport and thermoelectric 
properties of the AFM|| ground state as well as of the 
quasi-stable states AFMT, FMT, and FM||. Though the 
latter states correspond to energy larger than the ground- 
state one, they can be stabilized by some external forces, 
like a weak magnetic field, exchange coupling to an ad¬ 
jacent magnetic subsystems, etc. 


B. Phonon heat conductance 

To calculate some of the thermoelectric parameters, 
like thermoelectric efficiency (figure of merit), one needs 
to know not only the electronic contribution k to the 
heat conductance, but also the corresponding phonon 
term Kph. In a recent paper Yang et have deter¬ 
mined the phonon contribution to the heat conductance 
of narrow germanene and silicene nanoribbons, utiliz¬ 
ing the Landauer formula (ballistic transport approxi¬ 
mation) and the phonon transmission functions calcu¬ 
lated by an ab-initio method. We checked that /tph pre¬ 
sented in Ref.jl^ scales linearly with the nanoribbons 
width, and then extrapolated Kph up to zigzag nanorib¬ 
bons containing N = 8 zigzag chains (16 atoms in the 
zigzag chain), which are considered in this paper. We 
obtained the following values for N = 8 and T = 100 K: 
Kph = 0.92 nW/K for germanene, and Kph = 1.09 nW/K 
for silicene. Thus, for the model calculations we take 
Kph ~ 1 nW/K in Eqs (llOal) and (llObI) . 


C. Antiferromagnetic phase 

Now, we analyze ballistic transport and thermoelec¬ 
tric properties in the AFM configuration in more details. 
The corresponding electronic states for the AMFT (left 
panel) and AMF|| (right panel) configurations are pre¬ 
sented in Fig. [2] for three different values of the spin- 
orbit parameter. The corresponding transmission func¬ 
tions are presented in Fig. [31 The lowest value of Ago 
assumed in these figures is comparable to the spin-orbit 
parameter typical for silicene^ (note the factor 3-\/3 in 
the definition of spin-orbit coupling in this reference). 
However, larger values of the parameter Ago can appear in 
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FIG. 2: Electronic structure of the AMF_L (left panel) and 
AMF|| (right panel) states for indicated values of the spin- 
orbit parameter. Vertical dotted lines indicate the interval in 
the wavevector space, where the edge parameter ^ obeys the 
condition, |^| > 0.5. The other parameters are: U = 1.4 eV, 
V = 8, and t — 1.6 eV. 


other related two-dimensional crystals, like germanene^ 
(Ago = 43 meV) or stanene^ (Ago = 29.9 meV). More¬ 
over, results for larger values of the spin-orbit param¬ 
eter are very interesting from the point of view of the 
topological properties. Therefore, in numerical calcula¬ 
tions of the band structure, transmission, and thermo¬ 
electric coefficients we assume the following three values 
of Ago: Ago = 0.02t = 32 meV (strong spin-orbit cou¬ 
pling), Ago = O.Olt = 16 meV, and Ago = 1 meV (weak 
spin-orbit coupling). In turn, for the electron-electron 
interaction we assume [/ = 1.4 eV. 

To emphasize behavior of the edge states, which are 
responsible for topological properties, we displayed in 
Figs. [2] and |3] only a small energy region near the Fermi 
level Ep oi the corresponding charge-neutral nanoribbon. 
Moreover, only the region of wavevectors fc, in which the 
states are strongly localized near the left (L) or right 
(R ) edges of the nanoribbons is presented there. As 
follows from these figures, the electronic spectra reveal 
an insulating gap for small values of Ago, and this gap 
appears in both AFMT and AFM|| phases. Moreover, 


width of the gap is practically the same in both cases. 
Note, such a gap is absent in the limit oi U = 0, when 
zero-energy topologically protected states appear 4^ Fur¬ 
thermore, the valence and conduction states localized at 
the left and right edges are degenerate in both phases. 
These states are additionally spin degenerate. 
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FIG. 3: Transmission function for the AFMT (top panel) 
and AFM|| (bottom panel) phases, calculated for the same 
parameters as in Fig. 0 
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Qualitatively different behavior can be observed for 
larger values of Ago. The edge states are then still degen¬ 
erate, but in the AFMT phase they cross the Fermi level, 
so the system acquires properties of a topological insu¬ 
lator, whereas the AFM|| phase still exhibits properties 
of a conventional insulator. Thus, at a certain value of 
Ago, there is a transition from the conventional to topo¬ 
logical insulator in the AFMT phase, while no such a 
transition appears in the AFM|| phase. Moreover, one 
may conceive that by applying an appropriate external 
magnetic field it would be possible to switch between the 
AFM|| and AFMT states, and thus also to change the 
character of the system from conventional to topological 
insulator, and vice versa. Opening of an energy gap near 
the Fermi level at the transition from the AFMT state 
to the AFM|| can be clearly seen in the corresponding 
transmission function T{E) presented in Fig. [31 

Numerical results for the whole range of the angle a 
(not presented), i.e. including also the hypothetical un¬ 
stable magnetic states, show that for a fixed value of 
Ago the width of energy gap strongly depends on the an¬ 
gle a as well as on the Hubbard parameter U. Gener¬ 
ally, the gap increases with U and a, independently of 
























































































































6 



(i-Ep [meV] 


FIG. 4: Electrical conductance G (a) and electronic contribu¬ 
tion to the thermal conductance, ft, (b) in the AFM|| config¬ 
uration, calculated as a function of the chemical potential for 
indicated values of Aso and T = 100 K. The other parameters 
as in Fig.2. The insets show the corresponding results for the 
AFMT state. 


the nanoribbon width. Interestingly, the system is non¬ 
magnetic for a = 0 and U <1, but for higher values of a 
it exhibits the antiferromagnetic arrangement with quite 
considerable edge moments. Transition to the AFM state 
opens an energy gap, and the gap width monotonically 
increases with a. 

Thermoelectric properties in the AFMT phase were 
already analyzed in our earlier paper for a small value 
of Aso*^ Therefore, we focus here mainly on the AFM|| 
ground state, though some new features in the therm- 
electric properties in the AFMT state appear for large 
spin-orbit coupling, especially those associated with the 
transition from the conventional to topological insulator 
state. The energy gap in the AFM|| phase has a signifi¬ 
cant influence on the transport and thermoelectric prop¬ 
erties of the system. The electrical conductance G and 
the thermal conductance n due to electrons are practi¬ 
cally equal to zero in a wide energy region (correspond¬ 
ing to the gap) near the Fermi level Ep, and then they 
strongly increase at the gap edges, see Fig. S) This figure 
shows both G and fc as a function of the chemical poten¬ 
tial, and for both AFM|1 and AFMT (insets) phases. In¬ 
terestingly, both G and k in the AFM|| state only weakly 
depend on Ago. This is rather obvious as the correspond¬ 



g-Ep [meV] 


FIG. 5: Thermopower S for the AFM|| configuration, calcu¬ 
lated as a function of the chemical potential for indicated 
values of Aso and T = 100 K. The other parameters as 
in Fig.2. The inset shows the corresponding results for the 
AFMT state. 



FIG. 6: Figures of merit ZT for the AFM|| configuration, 
calculated as a function of the chemical potential for indi¬ 
cated values of Aso and T = 100 K. The other parameters as 
in Fig.2. The inset shows the corresponding results for the 
AFMT state. 

ing gap in the spectrum (see Figs. [5] and [3]) only weakly 
depends on the spin-orbit coupling parameter. Contrary, 
owing to the transition to the topological insulator state 
in the AFMT phase, and closure of the gap for large 
enough spin-orbit coupling, both G and k are constant 
and nonzero for chemical potentials in the region around 
Ep (see the insets in Fig. 01 It is also worth to note, that 
in some narrow regions of the chemical potential, the 
electrical and thermal conductances in the AFM|1 and 
AFMT phases are smaller for higher values of Ago- 
Due to the difference in the edge states, both AFMT 
and AFM 11 phases display also remarkably different ther¬ 
moelectric properties, especially the Seebeck coefficient 
(thermopower), see Fig. [5] Due to the gap in the AFM|| 
phase, see Fig. [31 the thermopower S considerably in¬ 
creases for chemical potentials in the vicinity of the Fermi 
level Fp, achieving very high values at low temperatures. 


























































7 


Since the gap weakly depends on the spin-orbit coupling, 
the thermopower also weakly depends on Ago- The situa¬ 
tion in the AFM_L state is different. Now, the gap disap¬ 
pears for large values of Ago, so the thermopower remains 
zero in a broad region of chemical potentials around 
see the inset in Fig. [S] Then, smaller peaks appear due 
to the narrow peaks in transmission around E k, ±200 
meV, as shown in Fig. [3] 

The presented results show that nanoribbons with rel¬ 
atively strong spin-orbit coupling and exhibiting AFM|| 
magnetic configuration may be attractive for applica¬ 
tions for conversion of thermal into electrical energy. The 
corresponding thermoelectric efficiency, described by the 
figure of merit ZT, is shown in Fig. |6l Unfortunately, 
ZT < 1, so the thermoelectric efficiency is below that 
anticipated for practical applications. However, one may 
conceive some ways of reducing the phonon heat conduc¬ 
tance in order to increase ZT. 


D. Ferromagnetic phase 

Now, we analyze ferromagnetic arrangement of the 
edge moments and compare the results with those ob¬ 
tained for the corresponding antiferroinagnetic phase, 
presented and discussed in the preceding section. Both 
FM± and FM|| configurations are quasi-stable states 
with the corresponding energies above the energy of the 
AFM|| ground state. The ferromagnetic configurations, 
however, can be stabilized externally by a magnetic field 
or due to proximity effects. 

The calculated band structure near the Fermi level Ep 
is presented in Fig. [7] for both FM± and FM|| phases, 
whereas the corresponding spin-dependent transmission 
functions are shown in Fig. [S] As in the AFM configura¬ 
tions, electronic spectra of both FM± and FM|| phases 
are similar for small values of Ago, and both phases are in 
the topological insulator state. The situation changes for 
larger values of Ago. According to Fig. [3 the FM± phase 
remains in the topological insulator state, independently 
of the spin-orbit coupling strength assumed in this fig¬ 
ure. However, for higher values of Ago the degeneracy of 
states localized at right and left edges of the nanoribbon 
is lifted. Significantly different qualitative behavior is ob¬ 
tained for the FM|| phase. Though the system is in the 
topological insulator phase for very small values of Ago, 
its character changes for stronger spin-orbit interaction, 
and the system transforms into a conventional insulator 
with a relatively wide energy gap. It is interesting to 
note, that the range of fc-wavevectors, where the edge 
states appear in the FM|| phase, is much narrower than 
in the corresponding FM± state as well as in the previ¬ 
ously discussed AFM phases (see the dotted vertical lines 
in Fig. [21). Moreover, the states localized at the left and 
right edges are roughly degenerate. 

The opening of an energy gap in the FM|| phase with 
increasing Ago and absence of such a gap in the FM± 
state are well seen in Fig. [51 where the spin-dependent 



K 71 K' K 71 K' 


FIG. 7: Band structure for the FM± (left panel) and FM|| 
(right panel) configurations. Vertical lines denote the interval 
of wavevector k, where the edge parameter ^ obeys the con¬ 
dition 1^1 > 0.5. The assumed parameters are: U = 1.4 eV, 
N = 8, and t = 1.6 eV. 


transmission function is presented as a function of energy. 
It is worth to note that the transmission function corre¬ 
sponding to a particular spin orientation, up or down, is 
asymmetrical near the Fermi level, i.e. transmission in 
the spin-up channel exhibits a peak above Ep^ while the 
appropriate peak in the spin-down channel appears below 
the Fermi energy. This practically applies to both FM|| 
and FM± configurations. The width of the energy gap in 
the FM|| state increases with increasing Ago. However, no 
gap appears in the FM± phase. Thus, the situation is re¬ 
markably different from that in the AFM configurations 
discussed above. 

In the ferromagnetic phases, the two spin channels con¬ 
tribute differently to electronic transport. If the spin re¬ 
laxation time is long and therefore spin mixing is negli¬ 
gible, the spin-up and spin-down channels can be treated 
as independent. One can then distinguish contributions 
of individual spin channels to the electrical conductance 
(and also to the Seebeck coefficient). The spin resolved 
electrical conductance is presented in Fig. [3^ as a func¬ 
tion of chemical potential. The wide central region with 
the conductance close to zero can be clearly seen for the 
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FIG. 8: Transmission in the FM_L (top panel) and FM|| (bot¬ 
tom panel) configurations, calculated for the same parameters 
as in Fig. 7. 


FM|| phase. This region corresponds to the energy gap. 
The conductance strongly increases near the gap edges, 
where it also remarkably depends on the spin orientation. 
In turn, the conductance in the FMT phase is finite and 
constant in the whole region of chemical potential shown 
in Fig. (see the inset), except the two peaks corre¬ 
sponding to the enhanced transmission in narrow regions 
of chemical potential for the spin-up {fj, > Ep) and spin- 
down {fi < Ep) channels, which appear near E Ri ±200 
meV. 

The thermal conductance due to electrons is shown in 
Fig. [91 d for two values of the spin-orbit parameter. Inter¬ 
estingly, in some regions of chemical potential the ther¬ 
mal conductance is smaller for strong spin-orbit coupling, 
see Fig. |9b. This appears mainly because the energy gap 
in the FM|| state is wider for stronger spin-orbit cou¬ 
pling, which effectively suppresses n in a broader region 
of the chemical potential. This difference is especially 
pronounced in the FMT state (see the inset in Fig. [SId), 
where both peaks, i.e. the one below and another one 
above the Fermi level, are notably reduced for the larger 
value of the spin-obit parameter Ago. 

In addition to the conventional thermopower, one can 
now also define the spin thermopower. Both the charge 
Sc and spin Ss Seebeck coefficients are presented in 
Fig. [TUlfor two values of the parameter Ago. Due to the 
energy gap in the FM|| state, the conventional Seebeck 
coefficient Sc is considerably enhanced in the vicinity of 
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g-Ep [meV] 



g-Ep [meV] 


FIG. 9: (a) Spin resolved electrical conductance G in the FM|| 
configuration for Ago = 32 meV. (b) Electronic contribution 
to the heat conductance n in the FM|| configuration for indi¬ 
cated values of Ago. The insets show the corresponding con¬ 
ductances in the FMT state. The results are for T = 100 K 
and other parameters as in Fig.7. 


the Fermi level. The enhancement strongly depends on 
Ago and increases for higher values of the spin-orbit cou¬ 
pling. Similarly, the spin thermopower Sg also exhibits 
a very pronounced peak for a narrow region of chemi¬ 
cal potential near Ep, and then remains rather small for 
chemical potentials outside this region. The peak in Sg 
is very narrow, much narrower than the relevant energy 
gap. Furthermore, height of the peak strongly depends 
on the spin-orbit parameter. Remarkably different be¬ 
havior has been obtained for the FMT phase which ex¬ 
hibits typical topological insulator features. Now, the Sc 
and Sg Seebeck coefficients are rather independent of the 
chemical potential in the vicinity of the Fermi level, as 
shown in the insets to Fig. [ini). The small peaks which 
appear for higher values of \^ — Ep\ follow from the nar¬ 
row peaks in the corresponding transmission functions, 
see Fig. [S] 

The corresponding figures of merit for the conventional 
and spin thermoelectricity are shown in Figs. Illh and 
lllb . respectively. Significant values of the figures of merit 
occur only in certain regions of the chemical potential, 
where the corresponding thermopowers reach local max- 
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H-Ep [meV] 


FIG. 10: Conventional thermopower Sc (a) and spin ther¬ 
mopower Ss (b) for the FM|| configuration and indicated val¬ 
ues of Aso- The insets show the corresponding thermopowers 
in the FMT state. The results are for T = 100 K, and other 
parameters as in Fig. 7 


ima. Both of them, however, are relatively small, simi¬ 
larly as it was in the case of AFM configurations. How¬ 
ever, especially small is the spin figure of merit, even if 
the corresponding spin thermopower is relatively high. 
This is because the spin thermoelectric efficiency in the 
narrow region of chemical potential around Ep, where 
the relatively high peak appears in the spin thermopower, 
is negligible (compare Fig. ITUb and Fig. Illbl. This sup¬ 
pression of ZTs appears due to vanishing spin conduc¬ 
tance in this region of chemical potential (note, the trans¬ 
mission functions for both spin orientations are equal this 
region, see Fig. [8]). 


IV. SUMMARY AND CONCLUSIONS 

We have analyzed the influence of Coulomb and spin- 
orbit interactions in nanoribbons of 2D buckled hexago¬ 
nal crystals with zigzag edges on the edge magnetic mo¬ 
ments and on the topological and thermoelectric proper¬ 
ties. We have shown that the Coulomb interaction, taken 
in the form of Hubbard term in the mean field approx¬ 
imation, leads to formation of the edge magnetic mo¬ 
ments, which for a sufficiently strong spin-orbit coupling 
are oriented in the nanoribbon plane in the ground state. 



FIG. 11: Figure of merit for the conventional thermoelectric¬ 
ity, ZTc, (a) and for the spin thermoelectricity, ZTs, (b) in 
the FM|| configuration and for indicated values of the spin- 
orbit coupling. The insets show the corresponding figures of 
merit for the FMT state. The results are for T — 100 K, and 
other parameters as in Fig.7. 


and the moments at one edge are opposite to those at 
the other edge (AFM|| phase) - in agreement with Lado 
et al^ For small spin-orbit coupling, this configuration 
is almost degenerate with the state, in which magnetic 
moments are oriented perpendicularly to the nanorib¬ 
bon plane (AFMT phase). States with ferromagnetic 
arrangements of the moments (FM|| and FMT phases) 
are quasi-stable and correspond to higher energies. Such 
states can be stabilized externally, e.g. by a magnetic 
field and/or due to proximity to a ferromagnetic contact. 

As the AFM 11 phase remains in the topological insu¬ 
lator state when the pin-orbit coupling increases, the 
AFMT phase undergoes then a transition from the topo¬ 
logical to conventional insulator, i.e. an energy gap opens 
in the corresponding edge states. This gap has a signif¬ 
icant influence on transport and thermoelectric proper¬ 
ties. First of all, the opening of the gap is associated with 
a nonzero thermopower in the vicinity of the gap edges. 
The ferromagnetic quasi-stable states, FM|| and FMT, 
behave differently with increasing Ago, as well. This be¬ 
havior qualitatively resembles that of the antiferromag¬ 
netic phases, but now the gap in the FM|| state, which 
opens with increasing Ago, is spin dependent. As a re- 








































































suit, a nonzero spin thermopower may be observed, in 2012/04/A/ST3/00372. 
addition to the conventional one. 
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